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ABSTRACT 



The inverse Compton scattering (ICS) model can explain various pulse pro- 
file shapes and diversity of pulse profile evolution based on the mechanism that 
the radio emission is generated through inverse Compton scattering between sec- 
ondary relativistic particles and radio waves from polar gap avalanches. In this 
paper, we study the parameter space of ICS model for 15 pulsars, which share the 
common pulse profile evolution phenomena that the pulse profiles are narrower 
at higher observing frequencies. Two key parameters, the initial Lorentz factor 
and the energy loss factor of secondary particles are constrained using the least 
square fitting method, where we fit the theoretical curve of the "beam-frequency 
mapping" of the ICS model to the observed pulse widths at multiple frequencies. 
The uncertainty of the inclination and viewing angles are taken into account in 
the fitting process. It is found that the initial Lorentz factor is larger than 4000, 
and the energy loss factor is between 20 and 560. The Lorentz factor is consistent 
with the prediction of the inner vacuum gap model. Such high energy loss factors 
suggest significant energy loss for secondary particles at altitudes of a few tens 
to hundreds of kilometers. 

Subject headings: model: Inverse Compton scattering model; methods: numerical; 
pulsars: general; pulsars: individual PSR B0301-M9, PSR B0525-F21, PSR 
B0540+23, PSR B0823+26, PSR B0919+06, PSR B1039-19, PSR B1133+16, PSR 

B1702-19, PSR B1706-16, PSR B2020+28, PSR B2021+51, PSR B2045-16, PSR 
B2154-F40, PSR B2310-F42, PSR B2319-F60) 
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1. Introduction 

Since the discovery of pulsars, a wealth of observational data have been accumulated 
for radio pulsars and the morphological characteristics of pulsar profile have been widely 
investigated. The core-double-cone model (Rankin 1983a,b), a widely used empirical model, 
attributes the various kinds of morphologies to geometrical origins. Assuming that the 
emission beam consists of a hollow core component, a nesting inner cone, and an outer 
cone, the observed single, double, triple, quadruple, and five-component pulse profiles can 
be explained as geometrical effects that the line of sight sweeps across the emission beam 
from different locations. Despite the success of the empirical model, physical mechanisms of 
radio emission, however, remain open. Among the models of coherent curvature radiation 
(Ruderman & Sutherland 1975, hereafter RS75, Gil & Snakowski 1990), plasma instabilities 
(Asseo et al. 1990, Luo & Melrose 1995, Weatherall 1998, Gedahn et al. 2002, Melrose 
& Luo 2004), and inverse Compton scattering (hereafter ICS, Qiao & Lin 1998, Qiao et 
al. 2001), the ICS model manifests itself in simplicity to generate the frequency-dependent 
beam structures and flexibility to reproduce various kinds of frequency dependence of 
the profile shape and the pulse width. Beside pulse morphologies, the high brightness 
temperature and polarization properties of pulsars can be also explained by the ICS model 
(Qiao & Lin 1998, Zhang et al. 1999, Xu et al. 2000). 

In the ICS model, radio emission is generated by the inverse Compton scattering 
process between secondary rclativistic electrons/positrons and initial low frequency waves 
(with frequency i/q ~ 10^~^ Hz), which are produced by avalanches of the inner vacuum 
gap. The relation between the altitude and the Lorentz factor of secondary particles 
determines the beam pattern. For example, if the Lorentz factor of secondary particles 
decreases as they flow out along open field lines, the emission at a given frequency would 
emerge at three different altitudes, which corresponding to the hollow core, inner and outer 
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cone components. If the Lorentz factor keeps constant, it would form a beam pattern of 
only one hollow core and one cone (the inner cone). The ICS model predicts that, in the 
outer cone, the emission with a higher frequency is generated at a relatively lower altitude, 
while the radiation in the inner cone follows the opposite relation. Observationally, the 
anti-correlation between the pulse width and the frequency is found in many pulsars with 
conal components, although a small amount of pulsars show constant or increasing pulse 
widths at higher frequencies. The ICS model agrees with these observations, and attributes 
these two types of relation to the outer and inner conal components, respectively (Qiao et 
al. 2001). 

In order to constrain the parameter space of ICS model, some authors have already 
compared the predictions of the ICS model with the observational data. There are 
three basic parameters for the ICS model, the initial Lorentz factor 70, the frequency of 
initial radio wave vq, and the energy loss factor ^ for the secondary particles. Lee et al. 
(2009, hereafter LCW09) found that the ratio between the radiation altitudes at different 
altitudes is insensitive to inclination angles a for radio pulsars with large linear polarization 
position-angle (PA) swing rate. Assuming that the Lorentz factor decreases with altitude, 
the authors constrained the parameter space of 70 and the energy loss factor ^ for five radio 
pulsars by fitting the ratios of profile width at multi frequencies. Four out of five pulsars 
show clear decreasing pulse width with frequency, while the other one shows nearly constant 
pulse width. 

In this paper, we extend the above reverse-engineering test to a larger sample of pulsars 
and check the parameter space of the ICS model, where 15 pulsars with anti-correlations 
between the pulse width and the frequency is selected to check the radiation behavior of the 
conal beam in the ICS model. Wc calculate the "beam- frequency mapping" of ICS model 
in Section 2. The geometrical method to calculate the beam width from the pulse width is 
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presented in Section 3. Details about data reduction and related results are given in Section 
4. Conclusions and discussions are summarized in Section 5. 

2. The beam-frequency mapping in the ICS model 

In the ICS model, low-frequency waves are excited by the periodic breakdown of 
the inner vacuum gap. The waves are then inversely Compton scattered by the relativistic 
secondary particles produced in the pair cascade process in the gap. The scattered waves 
have the frequency: 

1/ = 2721/0(1-1; cos ^i/c), (1) 

where v is the velocity of the particle and 9i is the incident angle between the particle 
motion direction and the incoming photons, and Uq is assumed to be 10^ Hz. There are a 
few physical considerations for the value of uq. RS75 pointed out that the growth and the 
fluctuation time scale of the inner vacuum gap might be (30-40)to, where tq is the time 
needed for conversion of a gamma-ray photon to a pair in the gap, which has the order 
~(gap height)/c. The fluctuation time scale is estimated to be about 10 microseconds. But 
it is not conclusive yet even with the recent simulation (Timokhin 2010). Observationally, 
durations of a few microsecond was indeed observed in some pulsars (e.g. Bartel & Sieber 
1978, Lange et al. 1998). 

For the geometrical conflguration in Figure 1, when the radiation altitudes are far from 
the pulsar surface, we have (Qiao et al. 2001) 

2cos^- (i?/r) (1-3 cos^e) 
cos Oi = y / n ) ^ 2 

7(1 + 3 cos2 9) [1 - 2{R/r) cos 9 + {R/rf] 
where r is the distance between the scattering point and the center of the pulsar, R is 
pulsar radius, 9 is the polar angle between the radiation location and the magnetic axis. 
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For a dipole magnetic field, we have 



r = Re sin 




(3) 



where Re is the maximal radius of a given magnetic field line. The angular beam width 
of the radiation beam coming from the place with polar angle 6 is, 

3 sin 2^ 
l + 3cos26' 

where 6^ is the angle between the direction of magnetic field at the point Q and the 
magnetic axis (Figure 1). From the radiation altitude r and Lorentz factor 7, one can 
determine the outgoing radio wave frequency u and the angular beam width 6^ using the 
above equations. 

The energy of the secondary particles decrease as they flow along the field lines, it is 
usually assumed that the Lorentz factor follows 



where 70 is the initial Lorentz factor at the top of the gap and ^ is the energy loss factor. 
From Equations (IT]) and ([5]), one can see that uq and 70 are degenerated, i.e. different 
groups of uq and 70 can give identical results in the ICS model. Thus, in the subsequent 
test for the ICS model, we fix the value of z/q without losing generality. 

Using the above equations, the relation between the outgoing frequency z/ and the 
beam width 9^, i.e. the function z/ = z/(z/o, 70, ^, 6^^), can be calculated. We refer this u-O^ 
relation as the "beam-frequency mapping" . A characteristic beam-frequency mapping with 
the dropping-rising-dropping pattern is given in Figure 2. Compared with the empirical 
core- double- cone framework for pulsar radiation, the leftmost branch corresponds to the 
core component, the middle rising part corresponds to the inner cone and the rightmost 
part corresponds to the outer cone. In this paper, we fit such beam-frequency mapping to 



7 = 7o[l-e(r-i?)/i?e], 



(5) 
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the observational data to infer the parameters of ICS model, where the fitting method is 
given in next section. 



Method 



In our method , we first determine Of, from the pu lse width using conventional 



geometrical model ( IGil 



1984 



Lyne fc Manchester 



19881 ). and then fit the beam-frequency 



mapping to 0^ data at multiple frequencies to infer the parameters. Here the pulse width 
is measured using the Gaussian decomposition method (Wu et al. 1992, Kramer et al. 
1994a,b, LCW09). 



3.1. Measurement of pulse width 

The Gaussian decomposition method is widely used to measure the pulse width of a 
pulse profile, where one models the profile with multiple Gaussian components and extracts 
the pulse widt h using curve fitting techniques, i.e. one fits the pulse profile with the 



following form (jWu et al. 



1992 



Kramer 



1994 : 



Kramer et al 



19941 ) 



exp 



i=l 



2af 



(6) 



Here we denote the intensity of pulse profile I {(f)) as a function of longitudinal phase 0, 
which is modeled with n Gaussians. The Jj is the amplitude of the i-th Gaussian, 0j and cFi 
are the central phase and standard deviation of Gaussian component respectively. 

We use the 10% width throughout this paper, which is defined as the full pulse width 
2A0 between the leading phase of the leftmost component and the trailing phase of the 
rightmost component down to the 10% level of their maximal intensities. 



The number of Gaussians, the amplitude, the peak phase and the typical width 
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of each Gaussian are free parameters. Following LCW09, we use Levenberg-Marquardt 
method to perform the least square fitting, which is accepted only when i) nonparametric 
Kolmogorov-Smirnov test is passed when comparing with the distributions of residues in the 
on-pulse and off-pulse regions, and ii) the difference between the rms levels of two residues 
is close to zero, i.e. rj = |cron — o'ofr|/o"ofr — 0. The number of Gaussians is determined 
by minimizing rj. We run a small Monte-Carlo simulation to determine the final fitting 
parameters, where the fitting is repeated several times (about 10) with randomly generated 
initial values, and the final parameters are the averaged values from each individual 
accepted fitting. We list the measured final pulse width in Table 1. 



3.2. Determining angular sizes of radiation beams 



From the pulse width, we calculate the angular size of radiation beam following 
the conventional geometrical model, which assumes that, i) the magnetic field in the 
magnetosphere is dipolar, and ii) the radiation direction is parallel to the local magnetic 
field (Figure 3). When the line of sight locate s in the Q — u plane, the rate of linear 



polarization PA swing reaches the maximum k (iRadhakrishnan &: Cooke 



1969I ). where 



sm a 
sin (3 



(7) 



Here a is the inclination angle, the impact angle /3 is defined as the angle between the line 
of sight and the magnetic axis. 



Given the inclination angle a and k, we can determine (3. We can further sub stitute 



A(j) and a into the 



Lyne fc Manchester 



oUow ing geometrical relation to calculate the beam width 6^ (iGil 



1984 



19881) 



sm 



2 / 



sm I I sm a sm[a + p J + sm I — 
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3.3. Measurement for ICS model parameter 

The ICS model parameters are derived via fitting the 6'^ data at multiple frequencies 
with the beam-frequency mapping. On the one hand, given frequency bands at 
central frequencies of z/j, we can measure N beam widths ^^^j from the observation, where 
i = 1 ■ ■ ■ N . On the other hand, with Equations (l)-(4), we can calculation the predicted 
beam width 9^^^^ of the ICS model at those frequencies as a function of parameters 70 and ^. 
Thus we can fit the ICS model to the data by minimizing the following standard x^, 

N / n /Dies \ 2 

where S6^ is the error of 0^. In our fitting, the third ICS model parameter, the background 
wave frequency uq, is fixed as 10^ Hz. Such fixation is due to the parameter degeneracy, 
i.e. the effect of z/q is completely absorbed into the parameter 70 as indicated in the 
Equation (1). 



4. Data reduction and results 

The pulse profile data are from the European Pulsar Network Database, where 15 
pulsars are selected according to two criteria, viz. 1) high pulse profile quality at more 
than five frequencies ranging over one order of magnitude, (2) the absolute values of 
the maximum slope rate of linear polarization PA n is much larger than 1. The second 



requirement, i.e. k ^ 1, is to make sure that the ratio between angu 



beams at different frequencies are insensitive to the inclination angle (iLee et al. 



ar sizes o: 



" radiation 



20091 ). 



As explained in the previous section, given the inclination angle a, the maximal PA 
swing slope k, and the pulse width, we can calculate the angular sizes 9^ of the radiation 
beams at multiple frequencies. With the beam-frequency mapping, we fit for the ICS model 
parameters. 
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One caveat is that it is difficulty to measure the inchnation angle accurately. 
Two types of methods have been used to estimate a. The ffist type of methods are 



fitting the polarization PA data with the conyentio nal (iRadhakrishnan fc Cooke 



1969 



Lyne fc Manchester 



1988 



Everett fc Weisberg 



20011 ) or modified versions (Blaskiewicz et 



al. 1991) of the rotation vector model. The second type of methods are based on some 
statistical relations between the pulse period and the opening angle of emission beam 
(Rankin 1990, Kijak & Gil 1997). Due to the limitation of each method and data quality, 
in many cases, different authors got inconsistent values of inclination angle. Two major 
methods are used here to reduce the effect of uncertainty of the inclination angle. Firstly, 
we select pulsars with larger k. As proved by LCW09, although the absolute beam width 
relies on the inclination angle, the ratio between the beam widths at multiple frequencies 
is insensitive to the inclination angle. Such invariant ratio is already able to determine the 
ICS model parameters. Secondly, we perform the measurement of 6'^ and the fitting with 
the ICS model for all the possible a values, of which the range is presented in Table 2. By 
enumerate all the possible a, we ensure that the invalid parameter space is really not viable. 

We show the allowed regions in the 70-^ parameter space in Figure HJ where the 
contours are at 2-cr level. We also show the result from the method in LCW09. The results 
immediately show that fitting the absolute values of can tell us how the allowed parameter 
space depends on the inclination angle. However, it is not sensitive to constrain 70, because 
the outer-cone branches of beam-frequency curves with the same ^ but different 70 are 
compressed very much along the 6'^- dimension, as presented in Figure. O The constrained 
parameter ranges are also listed in Table 2. The general features are summarized as follows. 

(1) Combining with all the samples, the allowed parameter space are 70 > 4000 and 
20 < ^ < 560. Since the Lorentz factor of secondary particles is likely well below 10^ 
according to inner vacuum gap model, we set the cutoff of 70 to be 20000 in calculation. 
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No clear correlation is found between the constrained parameter values and observational 
quantities, e.g. the pulse width, the surface magnetic field, the pulsar age and the spin-down 
energy loss rate. 

(2) It shows a general trend that both the initial Lorentz factor 70 and the energy loss 
factor ^ increase as a decreases. This is because 9^ and the emission altitude r become 
smaller when a decreases, it requires a faster energy loss so that the Lorentz factor decreases 
to proper values to produce the radio emission at the observed frequencies. 
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5. Conclusions and Discussions 

We collected a sample of 15 pulsars which have the anti-correlation phenomenon 
between the pulse width and the observing frequency. Their pulse widths are measured from 
multi-frequency profiles by using the Gaussian decomposition techniques. Beam widths are 
calculated with the classical beam geometry model for possible inclination angles. Then 
they are fitted with the beam-frequency relation of the ICS model to constrain two free 
parameters, i.e. the initial Lorentz factor 70 and the energy loss factor ^. The fitting is 
performed for a group of possible a. It shows a trend that 70 and ^ could be larger for 
smaller inclination angles. The allowed parameter space are 70 > 4000 and 20 < < 560. 
The range of the initial Lorentz factor is generally consistent with the prediction of the 
inner vacuum gap model (Timokhin 2010). The constrained values of ^ suggest that the 
secondary particles need to lose a large fraction of their initial energy. 

In our calculations, we assume z/q = 10^ Hz without losing generality due to the 
parameter degeneracy, which also agrees with the physical consideration of the inner 
vacuum gap model. If we take vq = 10^ Hz, according to Equation ([1]), the resulted Lorentz 
factors would be about 3 times higher than the present results. 

Our results indicates a bit higher Lorentz factor 70 > 4000 than the traditional picture 
of the inner vacuum gap model (~ 1000 in RS75), but this agrees with the requirement to 
produce the radio luminosity of pulsars in the ICS model. Noting that the primary particles 
usually has a Lorentz factor of 10^ ~ 10'' due to radiation reaction (e.g. RS75), the present 
results imply that the multiplicity of secondary particles is about a few hundred. With such 
a multiplicity, only a small fraction of secondary particles being coherent (Qiao & Lin 1998) 
or even incoherent radiation in some cases (Zhang et al. 1999) is sufficient to account for 
the observed radio luminosity, because the ICS emission of a single particle is much more 



efficient than the curvature radiation. iHarding fc Muslimov 



(120111 ) found that a slightly 
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asymmetric distortion can significantly increase the accelerating electric field on one side 
of the polar cap and, combined with a smaller field line radius of curvature, would lead to 
larger pair multiplicity. This increase of the primary accelerating electric filed may also be 
the origins of high Lorentz factors of the secondary particles. 

The result of 20 < < 560 indicates an efficient energy lose for the secondaries. It 
has been suggested that the resonant inverse Compton scattering between secondaries and 
thermal photons from the pulsar surface can cause efficient energy loss at a certain surface 
temperature, but it is only effective within about one stellar radius above the polar cap 
surface (Zhang et al. 1997, Lyubarskii & Petrova 2000). The other emission mechanisms 
are much less efficient (Sturner 1995) for secondaries, hence can be neglected. Therefore, 
the mechanism for such efficient energy loss is still unclear. LCW09 revealed that the decay 
of secondary Lorentz factor varies significantly from pulsar to pulsar. For some pulsars (e.g. 
PSR B2016+28), of which the profile width nearly keeps constant at different frequencies, 
a small value of ^ is thus inferred. Such diversity of energy loss behaviors may due to the 
interplay between the residual electrical field and the radiation reaction at higher altitudes, 
of which detailed studies is still demanded. 
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Table 1. Pulsar parameters and their 10% pulse widths as functions of observing 
frequencies 



Name P (s) B {10^^ G) k i/ (Hz) 2A(f){°) 77 Reference 



B0301+19 1.388 1.361 -17 



B0525+21 3.745 12.4 31 



B0540+23 0.246 1.97 -9 



B0823+26 0.531 0.96 -4.5 
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Table 1 — Continued 



Name P (s) B {10^^ G) k i/(Hz) 2A0 (°) 77 Reference 
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Table 1 — Continued 



Name P (s) B {10^^ G) k i/(Hz) 2A0 (°) 77 Reference 



B1133+16 1.188 



B1702-19 0.299 



B1706-16 0.653 



2.13 12 



1.13 -14 



2.05 -9 



0.102 22.9±1.3 

0.408 14.7±0.6 

0.925 13.4±0.7 

1.408 14.0±1.5 

1.71 12.3±0.8 

2.25 11.1±0.4 

4.85 11.2±0.6 

10.45 9.1±0.8 

0.61 20.3±0.7 

0.925 18.2±0.9 

1.408 17.4±0.4 

1.642 17.2±0.4 

4.85 16.7±0.4 

0.408 14.8±1.0 

0.61 14.2±0.6 

0.925 14.6±0.8 

1.408 14.5±0.5 

4.75 10.8±0.5 



1.633 


KL99 


0.095 


GL98 


0.077 


GL98 


0.959 


GL98 


0.614 


HX97 


U.D40 


T^V TOR 


2.485 


HX97 


0.577 


HX97 


0.388 


GL98 


0.121 


GL98 


0.437 


GL98 


0.407 


GL98 


0.091 


KKW97 


0.272 


GL98 


0.251 


GL98 


0.300 


GL98 


0.065 


GL98 


0.103 


SGG95 
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Table 1 — Continued 



Name P (s) B {10^^ G) k i/(Hz) 2A0 (°) 77 Reference 



B2020+28 0.343 



B2021+51 0.529 



B2045-16 1.962 



B2154+40 1.525 



0.82 15 



1.29 4 



4.69 -30 



2.32 8 



0.41 


18, 


.5±0 


0.925 


18, 


.5±0 


1.408 


19, 


.7±0 


1.642 


19, 


.4±0 


4.00 


1 s 
io. 


.D±U 


10.45 


14, 


.9±0 


0.41 


25, 


.6±1 


0.925 


22, 


.3±0 


1.408 


22, 


.9±0 


4. (0 


1 7 

1 1 . 


Q 1 n 
.oztU 


10.45 


13, 


.0±0 


14.6 


12, 


.7±0 


0.606 


17, 


.8±0 


0.925 


17, 


.9±0 


1.408 


16, 


.8±0 


1.642 


16, 


.4±0 


4.85 


15, 


.7±0 


0.408 


32, 


.1±2, 



5 0.522 GL98 

6 2.365 GL98 
4 0.549 GL98 
6 1.893 GL98 
8 1.780 HX97 
,7 0.275 HX97 

4 0.489 GL98 

5 0.788 GL98 
8 1.567 GL98 

3 0.177 SGG95 

4 0.045 HX97 
,6 0.027 KXJ96 

2 2.131 GL98 
8 0.080 GL98 

3 1.746 GL98 

4 1.788 GL98 

5 0.293 HX97 

1 0.410 GL98 
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Table 1 — Continued 



Name 


P(s) 




K V (Hz) 


2A0(°) 


7] 


Reference 








0.61 


31.4±0.7 


0.294 


GL98 








0.925 


28.8±1.4 


0.105 


GL98 








1.408 


29.9±2.0 


0.108 


GL98 








1.642 


27.4±1.2 


0.191 


GL98 








4.85 


26.3±1.2 


0.178 


HX97 




U.349 


U.2U 


A f\ A AO 

4 0.408 


19.4±1.2 


0.049 


T no 

GL98 








0.61 


18.4±0.9 


1.417 


GL98 








0.925 


19.6±1.1 


0.194 


GL98 








1.408 


17.8±0.7 


1.581 


GL98 








1.615 


17.6±0.5 


1.419 


SGG95 








4.85 


15.7±0.3 


0.221 


HX97 








10.55 


15.2±0.8 


0.127 


SGG95 


B2319+60 


2.256 


4.03 


-8 0.925 


24.7±0.8 


0.079 


GL98 








1.408 


23.6±1.0 


0.379 


GL98 








1.642 


22.6±0.3 


0.188 


GL98 








4.85 


19.8±0.6 


0.146 


HX97 








10.55 


20.6±1.1 


0.554 


SGG95 
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References. — Pulsar parameters of t he periods and the surfa ce magnetic field intensities 



are taken from A' 
curves k are from 



NF pulsar catalogue ([Manchester et al. 



19931 ) ■ The maximal slops of PA 



Lyne fc Manchesten (119881 ) . The original references for the profile data are 



given in the last column of the table as the following abbreviations: ANT94-Arzoumanian 
et al. 1994, GL98-Gould & Lyne 1998, HX97-Hoensbroech & Xilouris 1997, KKW97-Kijak 
et al. 1997, KL99-Kuz'min & Losovskii 1999, KXJ96-Kramer et al. 1996, SGG95-Seiradakis 
et al. 1995. 



-23- 



Table 2. Range of the inclination angle a and the constrained joint parameter space. 



Name a{°) Range of ^ Range of 70 References for possible a. 



B0301+19 


5-90 


50-420 


5000-20000 


1 


2 3 5 


B0525+21 


14-67 


50-560 


5500-20000 


1 


4 5 6 


B0540+23 


17-90 


10-90 


4000-16000 


1 


2 4 8 


B0823+26 


60-90 


26-50 


5800-10000 


1 


2 5 6 


B0919+06 


58-90 


34-80 


4200-7600 


1 


8 


B1039-19 


23-31 


175-410 


11000-20000 


5 


10 


B1133+16 


33-60 


80-425 


6500-20000 


1 


5 6 8 


B1702-19 


75 


64-110 


10000-20000 


5 




B1706-16 


40-70 


70-290 


6000-20000 


5 


6 


B2020+28 


60-71 


60-130 


9200-20000 


5 


6 


B2021+51 


16-30 


68-170 


11500-20000 


9 




B2045-16 


24-90 


70-375 


8000-20000 


2 


5 


B2154+40 


15 


210-430 


12000-20000 


5 




B2310+42 


24-56 


25-175 


4000-15000 


5 


6 7 


B2319+60 


12-35 


85-245 


10800-20000 


5 


8 



References. — 1-Blaskiewicz et al (1991); 2-Hoensbroech & Xilouris (1997); 
3-Narayan & Vivekanand (1982); 4-Gould (1994); 5-Kuz'min & Wu (1992); 
6-Lyne & Manchester (1988); 7- Wang & Wu (2003); 8-Kijak & Gil (1997) ;9- 
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Rankin (1993); 10-Everett & Weisberg (2001) 
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Fig. 1. — Geometry of the inverse Compton scattering model (Qiao & Lin 1998). 
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Fig. 2. — The "beam-frequency mapping" in the ICS model. For the curves from bottom to 
top, the corresponding 70 are 4000, 6000, 8000, 10000, 15000, while ^ is fixed at 50. It shows 
that Ofj, at a given frequency in the outer-cone branch becomes insensitive to 70 when 70 is 
large enough. 
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Fig. 3. — The geometry of pulsar radiation beam. The hne of sight sweeps the beam from 
A to B and the radiation is observed at C (Lyne & Manchester 1988). 



- 28 - 



Fig. 4. — The joint parameter space of the initial Lorentz factor 70 and the energy loss 
factor C,- The grey areas represent the allowed parameter space derived at 95% confidence 
level through fitting the index of 6^ — u relation, which is insensitive to inclination angles 
(LCW09). The solid-line and dotted-line contours represent the 2a confidence regions of 70 
and ^ constrained by fitting the absolute values of 6^ for different inclination angles, where 
the symbol "+" stands for the best fit 70 and ^. For each pulsar, some inclination angles 
within the possible range given in literature are selected as examples, their values are shown 
within the contour regions. For a particular inclination angle, the joint parameter space is 
the overlapped region of grey area and the contour corresponding to the inclination angle. 
See Table 2 for details of the joint parameter space. 
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